Scaling Laws and Topological Exponents in Voronoi Tessellations 
of Intermittent Point Distributions 
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Voronoi tessellations of scale-invariant fractal sets are characterized by topological and metrical 
j^jfj' properties that are significantly different from those of natural cellular structures. As an example we 

r-{ i analyze Voronoi diagrams of intermittent particle distributions generated by a directed percolation 

process in 2 + 1 dimensions. We observe that the average area of a cell increases much faster with 
the number of its neighbours than in natural cellular structures where Lewis' law predicts a linear 
behaviour. We propose and numerically verify a universal scaling law that relates shape and size of 
the cells in scale-invariant tessellations. An exponent, related to the topological properties of the 
(— | \ tessellation, is introduced and estimated numerically. 

■ 

(D ■ PACS numbers: 68.90.+g, 64.60.Ak, 82.70.Rr 
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Planar random mosaics are often encountered in nature, as, for example, in cuts of biological tissues jl],|| and 
two-dimensional soap froths . Thorough investigations of the topological and metrical properties of a huge variety 
of natural mosaics show a surprising similarity despite of the fact that the molding forces are completely different ||. 
\ The Voronoi construction || allows to generate a mosaic from a set of arbitrarily distributed seeds. It assigns to 
3 ■ each seed a cell which is defined as the set of all points of the plane which are at least as close to this seed as to any 
other seed. Since the resulting Voronoi tessellation is a space-filling cellular structure it allows to introduce the notion 

1 . of neighbourhood, i.e., two seeds are neighboured if their cells share at least one side. Thus the Voronoi construction 
enables us to investigate the neighbour statistics by analyzing the topological properties of the corresponding tessel- 
lation. This approach has been successfully applied to experimental monosize disc assemblies on an air table [Q,^), 
confirming the universality of random mosaics. In addition, the Voronoi construction offers an alternative tool for the 
characterization of the order-disorder transition which occurs when the density is reduced. 

■ So far all experimentally investigated cellular structures are characterized by a typical scale, namely the mean size 



of a cell. The question therefore arises how the mentioned universal laws are affected in the case of scale-invariant 
mosaics. In scale-invariant tessellations one expects that the area distribution of the cells exhibits a power-law 
Q\ ' behaviour. Moreover, it was shown in Ref. |l that in Sierpinski cellular structures the distribution of edge numbers 
t-H , also behaves algebraically. In the present work we are particularly interested in the relation between topological and 
metrical properties of scale-invariant mosaics. For this purpose we investigate Voronoi tessellations of intermittent 
particle distributions with algebraic long-range correlations that are generated by planar stochastic point processes. 

As an example we consider a Directed Percolation (DP) process in 2 + 1 dimensions Jl^,|ll| . In DP - interpreted as a 
time-dependent stochastic process - particles on a lattice either produce an offspring or self-destruct. Depending on the 
rates of offspring production and self-destruction this process exhibits a phase transition from a fluctuating active phase 
into an absorbing state without particles from where the system cannot escape. We use directed bond percolation |Tc| ] 
which is controlled by a single parameter, namely the percolation probability p. Below a critical threshold p < p c 
the system approaches the absorbing state in exponentially short time whereas in the active phase with p > p c a 
fluctuating stationary state exists on the infinite lattice. Close to criticality such a stationary DP process evolves 
through configurations that are characterized by spatially intermittent patterns with long-range correlations |l2| ], i.e., 
very few particles form highly localized clouds separated by large empty regions. 

In the present work we use spatial configurations (snapshots) of active particles generated by an almost critical 
stationary DP process as point sets for a Voronoi construction and to study the properties of the resulting tessellations. 
As shown in Fig. [|, such a tessellation appears to be very different from natural cellular structures. The reason is 
that point sets generated by an almost critical DP process approximate a fractal set, apart from lower and upper 
cutoffs due to system size and lattice spacing. In contrast to natural cellular structures the corresponding Voronoi 
tessellations are invariant under rescaling. Here we address the question how the topological and metrical properties 
of scale-invariant tessellations differ from those of natural structures. To this end we focus on the correlation between 
shape and size of the cells. 

The simplest quantity describing these correlations is (A)k/(A), the normalized average area of an arbitrarily 
chosen fc-sided cell. As pointed out first by Lewis pL {A)k increases linearly with k for a huge variety of tessellations. 
Deviations from Lewis' law have so far been observed in soap froth experiments |||J and air table tessellations in 
the dense packing regime 0. These deviations, however, are restricted to cells with only few neighbours whereas for 
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large fc the linear law appears to be asymptotically valid |4] 




FIG. 1. Ordinary and scale-invariant cellular structures. The figure shows parts of Voronoi tessellations for a Poissonian 
distribution of points (left) and for a point set generated by an almost critical DP process (right). 

In the present work the distribution P(k, A) of fc-sided cells with area A generated by an almost critical DP process 
is studied numerically. To this end a directed bond percolation process with parallel updates is simulated on a 
500 x 500 square lattice with periodic boundary conditions. The simulations are performed in the active phase with 
a small reduced percolation probability p — p c = 10 -4 , using the critical percolation threshold p c — 0.287338 |L3| |. 
After equilibration over 10 5 time steps the system evolves in a practically stationary state with on average about 
10 3 particles. The value of p — p c and the large lattice size ensure that the systems stays in the active phase. At 
intervals of 100 time steps the actual configuration of active particles is converted into the corresponding Voronoi 
tessellation 0]. Averaging over 1000 of such tessellations we determine the relative frequency P(k, A) as well as the 
mean area (A)k of a fc-sided cell. 
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FIG. 2. (a) Normalized average area (A)k/{A) of fc-sided cells generated by a DP process in comparison with a Poisson 
distribution [15], a cucumber [1], and a soap froth [3]. The inset shows a magnification for small values, (b) First three moments 
((fc — ko) q )A of the edge number distribution for cells with area A. The shift fco = 3 is explained in the text. 
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FIG. 3. Numerical results from Monte-Carlo simulations in the scaling regime, (a) Double-logarithmic representation of 
the relative frequency P(k, A) of cells with k edges and area A in units of squared lattice spacing. The dashed line indicates 
the location of the maxima, (b) Collapse of the same data using the scaling form ([^visualizing the scaling function $. The 
best collapse is obtained for the topological exponent r = 0.20(4). 

As shown in Fig. |]a, the curve for (A)k/(A) shows an untypical behaviour as it increases significantly faster than 
in cellular structures with a typical scale p^Jl|,p|. It neither resembles Lewis' linear law nor a quadratic behaviour 
that has been discussed in the context of a perimeter law for metallurgical grain structures Q] . The numerical results 
instead suggest an algebraic increase. In order to check the quality of the scaling, we measured the first three moments 
of the edge number distribution of cells with area A (see Fig. ||b). We observed that the best results are obtained if 
the edge numbers are measured with respect to their natural minimum fco = 3. The slopes are estimated by 0.13(2), 
0.26(3), and 0.38(5) for q — 1,2,3, respectively, indicating simple scaling. 

Our measurements for P{k, A) are shown in Fig. ||a. As A increases, the maximum of the curves is shifted along 
the dashed line to higher edge numbers k. The observation that - despite of the discreteness of k - all curves roughly 
have the same shape supports the conjecture that the cell topology in DP may exhibit universal scaling. 

Let us recall the scaling properties of DP. As usual, we denote the directed dimension (time) by the index 1 1 and 
the other spatial dimensions by _L. In the active phase close to criticality the stationary particle intensity p and the 
spatial and temporal correlation lengths £j_ , £i scale as |t(| 



P~(p-Pc) f3 , €±~{p-Pc)~ 



(P-Pc)-^ 



(1) 



The exponents /3, z/j_> aud v\\ are the three basic critical exponents of DP which usually determine other DP exponents 
by simple scaling relations. In 2 + 1 dimensions they have been numerically estimated [jllj by (3 — 0.584, v±_ — 0.734, 
and = 1.295. 

The Voronoi construction maps each active site onto a cell with a certain area A and a topology k. The cell topology, 
i.e., the number of neighbours, can be regarded as an additional quantity for the description of a DP process in more 
than one spatial dimension. In order to find an appropriate scaling relation, we conjecture - in analogy to Eq. (Q) 
- that the cell topology exhibits scaling properties similar to those of distances and densities. This means that the 
corresponding scaling regime is restricted to k < k, where k diverges close to the transition as 



K ~ (p-Pc)~ 



(2) 



Here n limits the scaling regime of edge numbers similarly as £j_ limits the scaling regime of distances in the active 
phase, r is a critical exponent which will be referred to as topological exponent of DP. At present it is not clear 
whether r is a novel exponent or related to the bulk exponents /3, by a scaling relation (as, for example, the 

critical initial slip exponent of DP pi). 

In the scaling regime we expect the probability distribution P{k, A) for cells with topology k and area A to obey 
the scaling form 



P(k,A) - A-' 11 <S>(kA- 



12) 



(k < k,A<£ 2 a 



(3) 
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where $(z) is a universal scaling function. The exponents 71,72 can be derived as follows: Summing over the edge 
numbers k we obtain the area distribution 



/>oo 

P(i)=rP(M)^ dkP(k,A) ~ A 72-71 . (4) 
ft -/o 

Assuming that < 72 — 71 + 2 < 1 we can compute the average area of cells for a given scaling length £j_ by 

(A) ~ [ X diAP(A) ~ ; (5 ) 

where a is the lattice spacing and d — 2 the spatial dimension. Since (^4) is inversely proportional to the intensity of 
points p ~ £j_ Wi we obtain the relation 71 = 72 + 2 — (3/dvi_ by comparing the exponents. On the other hand, the 
argument of the scaling function $ in Eq. (||) should be invariant under rescaling A — > b d A, k — > b T ^ u± k, leading to 
72 = r /dv±. Thus the expected scaling relation reads 



^^~v^~~ ~~v~~~ ^v—v-^w^ ^^v-^ 

P(fc, A) - ^r-fi)/du ± $(kA- T / du± ) , (fc < k, A < Ci) 



(6) 



In order to verify this scaling law, we numerically estimated the topological exponent t by data collapse (see Fig. ||b). 
Again, the best collapse is obtained if one replaces k by k — kg with fco = 3, which is the natural minimum for the 
number of sides of a cell. Our best estimate is r = 0.20(4), corresponding to 71 = 1.73(3) and 72 = 0.13(3). The 
collapse is fairly convincing, supporting the above scaling hypothesis. In addition, the result is in agreement with 
the measurements in Fig. 0b. We note that this result should not depend on the specific choice of the DP dynamics 
(e.g., bond or site percolation). Rather we expect t to be a universal exponent characterizing the DP class in 2 + 1 
dimensions. The question whether r is independent or related to the other scaling exponents is still open. 

It should be noted that the proposed scaling hypothesis needs further verification as the scaling regime of k extends 
over one decade only. The main limitation is the small value of r, i.e., each decade of scaling range in k requires 
about five decades of scaling range in A. However, it can be shown that Eq. (|J) holds exactly in the case of a simple 
Sierpinski gasket which strongly supports our scaling hypothesis |17j|. The scaling form assumes that the underlying 
point distribution is a simple fractal. This is the case for spatial cuts of critical DP states, as stated in Ref. jL2|. For 
multifractal cellular structures, however, the proposed scaling relation should be replaced by appropriate multiscaling 
laws 0. 

We thank W. Wolf for teaching us how to use the LEDA software package. 
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